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Abstract. The AGILE 7-ray satellite provides large sky exposure levels (>10^ cm^ s per year 
on the Galactic Plane) with sensitivity peaking at E ^IQQ MeV where the bulk of pulsar energy 
output is typically released. Its ~1 /is absolute time tagging capability makes it perfectly suited 
for the study of 7-ray pulsars. AGILE collected a large number of 7-ray photons from EGRET 
pulsars (>40,000 pulsed counts for Vela) in two years of observations unveiling new interesting 
features at sub-millisecond level in the pulsars' high-energy light-curves, 7-ray emission from pulsar 
glitches and Pulsar Wind Nebulae. AGILE detected about 20 nearby and energetic pulsars with good 
confidence through timing and/or spatial analysis. Among the newcomers we find pulsars with 
very high rotational energy losses, such as the remarkable PSRB 1509-58 with a magnetic field 
in excess of lO'^ Gauss, and PSR J2229H-61 14 providing a reliable identification for the previously 
unidentified EGRET source 3EG2227H-6122. Moreover, the powerful millisecond pulsar B1821- 
24, in the globular cluster M28, is detected during a fraction of the observations. 
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INTRODUCTION 

Poor 7-ray pulsar statistics has been a major difficulty in assessing the dominant mech- 
anism which channels pulsar rotational energy into high energy emission as well as 
understanding the sites where charged particles are accelerated. The large field of view 
of the AGILE Gamma-Ray Imaging Detector (GRID) [1] allows long uninterrupted ob- 
servations and simultaneous monitoring of tens of nearby radio pulsars belonging to the 
"7-ray pulsar region" of the P-P diagram characterized by 5 > 2 x 10^^ G and spin- 
down energy £rot > 1-3 x lO^'^ erg s ^ [2]. Here we present the results of three years of 
pulsar observations with AGILE. 

AGILE OBSERVATIONS AND TIMING 

Pulsar data were collected since the early phases of the mission. Timing observations 
suitable for pulsed signal analysis started in July 2007. AGILE pointings consisted of 
long exposures slightly drifting with respect to the starting pointing direction. At the 
end of October 2009 AGILE started observing in spinning mode due to reaction wheel 
failure, which is not affecting AGILE/GRID sensitivity and pulsar observations. Data 
screening, particle background filtering and event direction and energy reconstruction 
were performed by the AGILE Standard Analysis Pipeline. We adopted the AGILE event 



extraction criteria and timing procedures calibrated and optimized with the observations 
of known y-ray pulsars as described in [3]. 

Given the measured AG/L£"s time tagging accuracy of <200 jUs and the good radio 
monitoring (i.e. valid epoch range, adequate number of time of arrivals (ToAs)) available 
for the majority of our targets, the most significant pulsed signal detection is typically 
expected within the errors of the radio ephemeris frequency values. In particular, we 
performed standard epoch folding and hen adequate radio observations covering the time 
span of the y-ray observations are available (i.e. where the WAVE terms are included in 
TEMP02 ephemeris files), we account also for the pulsar timing noise in the folding 
procedure as reported in [3] over a frequency range defined according to 3 a errors of 
radio ephemeris. Pearson's x'^ statistics is applied to the 10-bin folded pulse profiles 
resulting from each set of spin parameters, yielding the probabilities (weighted for the 
number of trials performed on the data set) of sampling a uniform distribution, assessing 
the significance of the pulsed signal (sinusoidal or not). Furthermore, we verified our 
timing results also applying bin-independent parameter-free statistics as the test and 
the //-test that are typically more sensitive than tests for the search of sinusoidal 
pulses. Firstly, for each target we searched for pulsed signal using the whole available 
data span. Later, each observation block was analysed to check for possible flux and/or 
pulse profile variability. We also performed a maximum likelihood analysis (ALIKE 
task) on the AGILE data for the regions containing our targets in order to exploit the 
instrument's imaging capabilities to assess y-ray source parameters. Here, we will focus 
predominantly on the timing analysis. 

In order to perform AGILE timing calibration through accurate folding and phasing a 
dedicated radio pulsar monitoring campaign (that will continue during the whole AGILE 
mission) was undertaken, using the telescopes of the European Pulsar Timing Array, as 
well as the telescopes of Parkes in Australia and Mt Pleasant in Tasmania. Geminga is a 
radio-quiet pulsar whose ephemeris can be obtained from the regular X-ray observations 
by XMM-Newton. In order to verify the performances of the timing analysis procedure 
described above, a crucial parameter to check is the difference between pulsar rotation 
parameters derived from radio. X-ray and y-ray data. The implementation of the folding 
method described in [3] allowed for a perfect match between the best period resulting 
from y-ray data and the period predicted by the radio ephemeris with discrepancies 
■^Crab ~ 3 X 10~^^ s. Comparable to the period search resolution rcnb ~ 2 x 10~^^ 
s. Ignoring timing noise in the folding process would yield discrepancies (and light 
curve smearing) which are expected to grow when considering longer observing time 
span. Thus the contribution of timing noise should be considered both in high-resolution 
timing analysis and in searching for new y-ray pulsars. 

"OLD" AND NEW y-RAY PULSARS 

The resulting high-energy light-curves for Vela, Crab, Geminga and PSR B 1706-44 
pulsars are shown in Figure 1 together with a sample of newly discovered y-ray pulsars. 

The availability of radio observations bracketing the time span of the y-ray obser- 
vations (or of X-ray observations very close to the y-ray observations for Geminga) 
allowed us to also perform accurate phasing of multi- wavelength light-curves. We found 




FIGURE 1. (Left) Light-curves of the /iGRET pulsars obtained during the first year of AGILE scien- 
tific operations. The dashed fine corresponds to the the radio main peak, (see [3] for details). (Right) 
Background-subtracted £ > 100 MeV folded pulse profiles (except for the E > 5Q MeV Jl 357-6429 
light-curve) of a sample of new y-ray pulsars discovered by AG/LZi. For PSR J2229H-6114, the £ > 1 GeV 
histogram is also shown. The main radio peak (1.4 GHz) corresponds to phase (see [4] for details). 



that the phasing of the AGILE light-curves of the four pulsars (radio/X-rays/y-ray peaks 
phase separations) is consistent with EGRET measurements (see [3] for details) imply- 
ing no evidence of systematic errors in absolute timing with an upper limit terr < 1 ms. 
The resulting fluxes (Pulsed Counts/Exposure) are also consistent with those reported in 
the EGRET Catalogue. 

The Crab SuperAGILE light-curve (ref) was produced with the same folding method 
reported for the GRID. Inspection of Figure 1 shows that the X-ray peaks are aligned 
with the E > 100 MeV data within A0 ~ 400 jUs (a value obtained fitting the peaks with 
Gaussians) providing an additional test of the AGILE phasing accuracy. 

The effective time resolution of the light-curves will obviously improve with expo- 
sure time At oc T~p and a resolution <50 jUs was obtained after two years of AGILE 
observations of Vela. 

Although AG/L£' multi-wavelength phasing of the four y-ray pulsars is consistent with 
the results obtained by EGRET, the plots shown in Figure 1 allow us to start assessing 
new features in y-ray pulsar light-curves. Narrower and better resolved main peaks are 
revealed, together with previously unknown secondary features. In particular, a third 
peak is possibly detected at ~3.7<7 level in the Crab light-curve (P3 in Figure 1) and 
some interesting features seem present in the Vela light-curves (confirmed by recent 
Fermi observations [5]). 

In any case, the highly structured light-curves hint at a complex scenario for the sites 
of particle acceleration in the pulsar magnetospheres, implying different electric gaps 
with physical properties probably mostly related to their height above the neutron star 
surface. Alternatively, slight spatial oscillations of the gap locations on timescales < 1 
day could be invoked to explain the multiple contiguous peaks seen in the light-curves. 

In this perspective, the AGILE light-curves time resolution, currently limited only by 
the (continuously increasing) source counts statistics, will eventually yield a pulsar gaps 
map (e.g. parametrized with adjustable accelerating electric fields strength, and location 
in the magnetosphere) by coupling timing analysis and phase resolved spectral analysis. 



The long monitoring of Vela pulsar (~40,000 pulsed counts in two years of observa- 
tions) allowed us to detect possible y-ray emission from pulsar glitches [3]. 

Furthermore, the Vela pulsar wind nebula was recently firmly detected by AGILE 
constraining the particle population responsible for the GeV emission and establishing a 
class of y-ray emitters that could account for a fraction of the unidentified galactic 7-ray 
sources [6] (see last Section) 

The radio-aligned light-curves of a subset of newly discovered y-ray pulsars for 
which our timing analysis yielded a >4o detection are also plotted in Figure 1. We 
note that in all cases radio and y-ray timing results are compatible, with the highest 
significance frequency detected in y-rays within the errors of the radio ephemeris value, 
considering also the period search resolution. Furthermore, we verified that our analysis 
procedure (potentially affected by instrument-related systematic errors and biases in 
events extraction criteria) does not produce fake detections at a significance level above 
3 a when the radio-ephemeris are applied to randomly extracted AG/LE data. The most 
significant detection is PSR J2229-I-6114 for which AGILE detected pulsed emission 
(radio/y-ray periods discrepancy <10^^ s) and pinpointed the most likely position to 
~0.2 deg from the radio pulsar. Our detection provides a reliable identification for the 
previously unidentified EGRET source 3EG 2227-1-6122. The AGILE source position, 
pulsed flux, and photon index (~2.2) are consistent with the EGRET values. The y-ray 
light-curve of this pulsar (detected up to over 1 GeV), featuring just one prominent peak 
shifted ~ 180 deg in phase from the radio main peak, is shown in Figure 1. 

PSRJ1513-5908 (B1509-58) was detected by COMPTEL in the 1-10 MeV range, 
while EGRET reported only marginal evidence for a weak <4o source at ~ 1 deg from 
the radio position, with a pulsed emission upper limit of <58 x 10~^ ph cm~^ s~^ (see 
dedicated section). 

At variance with all the other targets, the millisencond pulsar J 1824-2452 in the 
Globular Cluster M28 was detected hy AGILE, with good significance (>4a) and perfect 
radio-y periods match, only in the time interval 54339-54344 MJD. The main radio 
peak at 1.4 GHz is coincident with the broad single peak seen in y-rays. Only marginal 
detection was obtained integrating other observations with comparable exposure or the 
whole data span. Noise fluctuations could possibly explain the apparent variability. 
Alternatively, although its y-ray efficiency and high-stability of spin parameters are 
compatible with rotation-powered emission, some additional mechanism disturbing the 
neutron star magnetosphere in the dense cluster environment could be invoked to explain 
the variable y-ray phenomenology of this peculiar pulsar. AGILE timing failures at 
sub-millisecond level in some observations (mimiking source variability) cannot be 
excluded. However, this seems unlikely, since we verified both timing accuracy and 
stability at ~ 200 /is level with Vela pulsar observations. Confirmation of this tantalizing 
result about physical variability will rest on longer monitoring campaigns. 

THE "SOFT" PSR B1509-58: A NEW CLASS OF y-RAY PULSARS? 

PSR B 1509-58 has a period P= 150 ms and a period derivative 1.53 x lO-^^sg-i. 
assuming the standard dipole vacuum model, the estimated spin-down age for this pulsar 
is 1570 years and its inferred surface magnetic field is one of the highest observed for an 




FIGURE 2. (Left) Phase-aligned 7-ray light-curves of PSR B 1509-58 with radio peak at phase 0. 
(Right) SED of PSR B 1509-58 (solid Hne) obtained from a fit of pulsed fluxes from soft to hard 7- 
rays. The circular points represent COMPTEL observations. The square points show AGILE pulsed flux 
at 30 < £■ < 100 MeV and 100 < £ < 500 MeV. The horizontal bar represents AGILE upper Umit above 
500 MeV. The red arrow represent AGILE upper limit above 500 MeV. The two green arrows represent 
Fermi upper limits [7]. COMPTEL data are from [8]. 



ordinary radio pulsar: B ■^3.1 x 10 G, as calculated at the pole; its rotational energy 
loss rate is £ 1.8 x 10^^ erg/s. 

The young age and the high rotational energy loss rate made this pulsar a promising 
target for the 7-ray satellites. In fact, the instruments on board of the CGRO observed 
its pulsation at low 7-ray energies, but it was not detected with high significance by 
EGRET. AGILE obtained the first detection of PSR B 1509-5 8 in the EGRET band [4] 
confirming the occurrence of a spectral break. We then presented the results of a ~ 2.5 yr 
monitoring campaign of PSR B 1509-58 with AG/LE, improving counts statistics, and 
therefore lightcurve characterization, with respect to earlier AG/LE observations [9]. 

AGILE devoted a large amount of observing time to the region of PSR B1509- 
58. We observed PSR B 1509-58 in three energy bands: 30-100 MeV, 100-500 MeV 
and above 500 MeV. We did not detect pulsed emission at a significance a > 2 for 
E > 500 MeV. The 7-ray lightcurves of PSR B 1509-58 for different energy bands are 
shown in Fig. 2. The AGILE E > 30 MeV lightcurve shows two peaks. The first peak is 
coincident in phase with COMPTEUs peak [8]. In its highest energy band (10-30 MeV) 
COMPTEL showed the indication of a second peak (even though the modulation had 
low significance, 2.1a). This second peak is coincident in phase with AG/L£"s second 
peak (Fig. 2). AGILE thus confirms the previously marginal detection of a second peak. 

Based on our exposure we derived the 7-ray flux from the number of pulsed counts. 
The pulsed fluxes in the three AGILE energy bands were Fy = 10(4) x 10^^ ph cm^^ s~^ 
in the 30-100 MeV band, Fy = 2.1(5) x 10"^ ph cm^^ in the 100-500 MeV band 
and a la upper limit < 8 x 10"^ ph cm"^ s"^ for E > 500 MeV. 

Fig. 2 shows the SED of PSR B 1509-58 based on AGILE'S and COMPTEUs observed 
fluxes. COMPTEL observations suggested a spectral break between 10 and 30 MeV. 



AGILE pulsed flux at energies E > 30 MeV confirms the presence of a soft spectral 
break at E < 100 MeV As shown in Fig. 2, we modeled the observed fluxes with a 
power-law plus cutoff fit: F(E) =kx E~" exp[—{E / Ec)^], with three free parameters: 
the normalization k, the spectral index a, the cutoff energy Ec and allowing j8 to assume 
values of 1 and 2. No acceptable x'^ values were obtained for /3 = 2, while for an j8 = 1 
we found = 3.2 for v = 2 degrees of freedom, corresponding to a null hypothesis 
probability of 0.05. The best values thus obtained for the parameters of the fit were: 
k = 1.0(2) X 10-4, a = 1.87(9), Ec = 81(20) MeV. 

The bulk of the spin-powered pulsar flux is usually emitted in the MeV-GeV energy 
band with spectral breaks at < 10 GeV (e.g. [10]. PSR B1509-58 has the softest 
spectrum observed among y-ray pulsars, with a sub-GeV cutoff at £ «i 80 MeV. 

When PSR B 1509-58 was detected in soft 7-rays but not significantly at £■ > 30 MeV, 
it was proposed that the mechanism responsible for this low-energy spectral break might 
be photon splitting [12]. The photon splitting [11] is an exotic third-order quantum 
electro-dynamics process expected when the magnetic field approaches or exceeds the 
critical value defined as Bcr = m^c^ / (eh) = 4.413 x 10^^ G. In very high magnetic fields 
the formation of pair cascades can be altered by the process of photon splitting: 7 — > 77. 

In the case of PSR B 1509-58 a polar cap model with photon splitting would be able 
to explain the soft 7-ray emission and the low energy spectral cutoff, now quantified 
by AGILE observations. Based on the observed cutoff, which is related to the photons' 
saturation escape energy, we can derive constraints on the magnetic field strength at 
emission, in the framework of photon splitting: 

e,^,^~ 0.077(5' sin 0fe5,o)~^/^ (1) 

where £esc is the photon saturation escape energy, B' ~ B/Bcr and 9kB,o is the angle 
between the photon momentum and the magnetic field vectors at the surface and is here 
assumed to be very small: 0^-5 < 0.57° [12]. Using the observed cutoff (E = 80 MeV) 
we find that B' > 0.3, which implies an emission altitude < 1.3i?iV5> which is the height 
where also pair production could ensue. This altitude of emission is in perfect agreement 
with the polar cap models Additionally, PSR B 1509-58 [8, 13] shows evidence of an 
aligned geometry, which could imply polar cap emission. 

The polar cap model as an emission mechanism is debated. From the theoretical point 
of view, the angular momentum is not conserved in polar cap emission (see [14] in 
this book). And a preferential explanation of the observed 7-ray lightcurves with high 
altitude cascades comes from the recent results by Fermi [10]. 

Alternatively, an interpretation of PSR B 1509-58 emission can be sought in the frame 
of the three dimensional outer gap model [15]. According to these estimates a magnetic 
inclination angle a ^ 60" and a viewing angle ^ 75'' are required to reproduce the 
observed lightcurve. Finally, using the simulations of Watters et al. [16], the observed 
lightcurve from AGILE is best reproduced if a ^ 35° and ^ a; 90°, in the framework of 
the two pole caustic model. 

The values of a and ^ required by the model in [15] are not in good agreement with 
the corresponding values obtained with radio measurements. In fact, Crawford et al. [13] 
observe that a must be < 60° at the 3a level. The prediction obtained by the simulations 
in [16] better agrees with the radio polarization observations. In fact, Crawford et al. 
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FIGURE 3. Gamma-ray high- and very-high energy spectral distribution (v/v) of the Vela X PWN. 
HESS data fit an IC process (scattering on CMBR) related to electron power law index 2.0 with a 
break at 67 TeV and a total energy content of 2.2 x 10"*^ erg [24]. AGILE data are compatible with 
IC emission from the additional electron component well reproducing the observed total radio spectrum 
(iitot — 4-x 10^** erg) assuming the same ^ 5 /iG field strength as required by the TeV spectral break. 

also propose that, if the restriction is imposed that C, > 70° [17], then a > 30° at the 
3a level. For these values, however, the Melatos model for the spin down of an oblique 
rotator predicts a braking index n > 2.86, slightly inconsistent with the observed value 
(n = 2.839(3)). Therefore, at present the geometry privileged by the state of the art 
measurements is best compatible with polar cap models. 

PULSARS AND THEIR ENVIRONMENTS 

The Vela supernova remnant (SNR) is the nearest SNR (~290 pc) containing a bright 
pulsar, PSR B0833-45. The detection of Vela X at TeV energies demonstrated that this 
source emits non-thermal radiation, in agreement with the hypothesis that it corresponds 
to the pulsar wind nebula, displaced to the south by the unequal pressure of the reverse 
shock from the SNR [18]. The multiwavelength spectrum of the center of Vela X can 
be modeled as synchrotron radiation from energetic electrons within the cocoon (radio 
and x-rays) and inverse Compton emission from the scattering (by the same electron 
population) of the cosmic microwave background radiation (CMBR), the Galactic far- 
infrared radiation (FIR) produced by reradiation of dust grains, and the local starlight 
[19, 20]. Alternatively, a hadronic model can be invoked for the y-ray emission from the 
Vela X cocoon, where the emission is the result of the decay of neutral pions produced 
in proton-proton collisions [21]. Observations in the MeV-GeV band (HE) are crucial to 
distinguish between leptonic and hadronic models as well as to identify specific particle 
populations and spectra. 

With the aim of performing a sensitive search for close faint sources excluding the 
bright emission from the Vela pulsar, we discarded the time intervals corresponding to 
the phase interval 0.05-0.65. The analysis of the resulting off-pulse images unveiled 
few y-ray sources, none of which coincides with the Vela pulsar. A maximum likelihood 
analysis, performed on the E > 100 MeV dataset within a region of 5° around the pulsar 



position, revealed two sources at better than 3(7 confidence, the brightest of which, AGL 
J0834-4539 (~5.9a significance, -264 counts, Fy = (35 ±7) x 10"^ ph cm^^ g-i 
E > 100 MeV), is positionally concident with HESS 10835^55, the TeV source that 
is identified with the Vela X nebula [24], and has a similar brightness profile to it. This 
implies that AGL J0834-4539 is associated with the pulsar's PWN. 

The AGILE spectral points are a factor ~2 below the previous EGRET upper limits 
[22] and well above the extrapolation of HESS vFy spectrum to lower energies. In the 
frame of leptonic models, the AGILE measurements are not consistent with a simple 
multiwavelength spectral energy distribution involving a single electron population. The 
AGILE spectral points are one order of magnitude above the fluxes expected from the 
electron population simultaneously fitting synchrotron x-ray emission (peaking at ~1 
keV) and inverse-Compton (IC) TeV emission [19]. Additional electron populations 
should be invoked to explain the observed GeV fluxes. This is not surprising in view of 
the complex morphology of the PWN seen in radio and x-rays, where different sites and 
features of non-thermal emission are present. In fact, the position where AGILE sees the 
maximum brigthness is also roughly where the 8.4 GHz radio emission is brightest [23]. 
AGILE data are compatible with the IC parameters modelled by de Jager et al. (electron 
spectral index 1.78 and maximum energy —20 GeV), although our measurements could 
suggest a higher contribution from IC photon seeds. AGILE measurements would be 
incompatible with the scenario of nucleonic y-ray production in the Vela TeV nebula in 
the frame of a single primary electron population. 
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